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We simulated the classical two-dimensional aniso- 
tropic Heisenberg model with full long range dipole 
interaction with an algorithm especially designed 
for long range models. The results show strong ev- 
idence for a first order reorientation transition at a 
temperature Tr < Tq for appropriate parameters 
of the model Hamiltonian. 

Introduction 

Recent experiments on thin ferromagnetic films 
have shown that the magnetic behaviour of these 
systems is strongly affected by anisotropy fields 
which act on the spins in the surface [1]. These 
experiments show that for many systems the direc- 
tion of the magnetization vector is a very sensitive 
function of temperature. At low temperatures the 
anisotropy fields force the magnetization vector to 
be perpendicular to the surface, but with increas- 
ing temperature the magnetization vector switches 
to an in plane direction. The switching tempera- 
ture decreases with increasing film thickness. This 
reorientation transition has been investigated the- 
oretically using classical [2] and quantum mechani- 
cal [3] meanfield theory and renormalization group 
theory [4]. These methods predict a reorientation 
transition of either first order in classical systems or 
second order in quantum systems, while a low tem- 
perature series expansion shows no evidence for a 
reorientation transition [5]. These conflicting re- 
sults call for further investigations. 

The model 

Therefore we have performed Monte-Carlo simu- 
lations on a classical two-dimensional anisotropic 
Heisenberg system with long range dipole interac- 
tions on a square lattice. The model is described 
by the Hamiltonian 
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where the s; = (s?,s?,sf) are classical Heisenberg 
spins at position r; = (x t ,y t ,0) and r t j = r t — 
Tj. J is the nearest-neighbour exchange coupling 
constant, K is the on-site anisotropy, and D is the 
strength of the long range dipole interaction. We 
consider ferromagnetic systems with J = 1 in this 
paper. 

The ground state energy per spin of a periodic sys- 
tem of fixed size can be calculated exactly as long 
as D is sufficiently small compared to J to ensure 
that all spins are aligned parallel. Using the two- 
dimensional lattice sum [6] 

oo 
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where a) is the generalized Riemann zeta func- 
tion and Re(s) > 1, one obtains the ground state 
energy per spin as a function of the z-component 
of the ground state magnetization Sq S as 

eosisbs) = -2J-K(s* GS f + ^ (3(^ s ) 2 -l)6(^). 

Hence all spins are pointing in z-direction when 

Ti < 30P72) = °- 1476 > while for Ti > 30P72) the 
spins are oriented in an arbitrary direction in the 
xy-p\a,ne. 

Method 

In the Monte-Carlo simulation we considered qua- 
dratic systems with periodic boundary conditions 
and summed the long range dipole interaction over 
all periodic images of the system. The simulations 
were done with an algorithm especially designed 
for long range systems: The local fields at each 
site are computed at the beginning of the simula- 
tion and are only updated when a spin flip attempt 
(SFA) is accepted. With this method we obtain a 
speedup factor of approx. 30 in the relevant tem- 
perature regime. To sample the phase space more 
efficiently, we combine both local and global spin 
update mechanisms: Single Heisenberg spins are 
chosen randomly and are reoriented into a random 
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Figure 1: Components of the total magnetization as 
a function of temperature for a 16 x 16 system with 
K = 0.2, D = 0.028 
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Figure 2: Probability density determined by the simu- 
lation /i(|s z |) of the absolute value of the z-component 
of the total magnetization at several temperatures for 
a 32 x 32 system with K = 0.2, D = 0.028 



direction. Additionally, after a fixed number of lo- 
cal SFAs, a global SFA is made, where all spins are 
mirrored at a plane with randomly chosen normal 
vector. Note that this transformation leaves the 
exchange part of the Hamiltonian invariant. This 
global spin flip is similar to the cluster spin flip 
in the Swendsen-Wang Algorithm for Heisenberg 
models [7], although this algorithm is not applica- 
ble in the present model due to the anisotropy of 
the dipole interaction. 



Results 

In figure 1 the perpendicular component (s z ), the 
in-plane components (s xy ) and the absolute value 
of the spontaneous magnetization for a system 
of 16 x 16 spins are depicted. The values are ob- 
tained from the maxima of the distribution function 
described below. In this run we made 512000/site 
local and 2000 global SFAs per temperature step 
using 60 CPU-hours on an IBM RS/6000 590 work- 
station. As one can see clearly, the magnetization 
changes from perpendicular to inplane orientation 
in the temperature range 0.12 < Tr < 0.13 for this 
set of parameters. 

To investigate the order of the reorientation tran- 
sition, we considered the probability density 
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h(x) = lim - y 6(x T — x) 

T = l 

for x = \s z \, as determined quite accurately from a 
histogram for a 32 x 32 system at fixed tempera- 
ture. Here a typical run involved ca. 2 x 10 6 local 
SFAs per site and 7500 global SFAs and took ss 10 
CPU-hours on an IBM 590. h(\s z \) is depicted in 
figure 2. One can see that the reorientation transi- 
tion seems to be of first order, as the distribution 
has two maxima at Tr. 



Conclusions 

To summarize, we simulated the classical two- 
dimensional anisotropic Heisenberg model with full 
long range dipole interaction with an algorithm es- 
pecially designed for long range models. The first 
results show strong evidence for a first order re- 
orientation transition at a temperature Tr < Tc- 
This phase transition will be investigated in further 
detail in the future. 
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